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Abstract. We consider the Gaudin model associated to a point z G C" with pairwise 
distinct coordinates and to the subspace of singular vectors of a given weight in the 
tensor product of irreducible finite-dimensional s?2-representations, |^. The Bethe 
equations of this model provide the critical point system of a remarkable rational sym- 
metric function. Any critical orbit determines a common eigenvector of the Gaudin 
hamiltonians called a Bethe vector. 

In [HeV] . it was shown that for generic z the Bethe vectors span the space of singular 
vectors, i.e. that the number of critical orbits is bounded from below by the dimension 
of the space of singular vectors. The upper bound by the same number is one of the 
main results of jSVj . 

In the present paper we get this upper bound in another, "less technical" , way. The 
crucial observation is that the symmetric function defining the Bethe equations can be 
interpreted as the generating function of the map sending a pair of complex polynomials 
into their Wronski determinant: the critical orbits determine the preimage of a given 
polynomial under this map. Within the framework of the Schubert calculus, the number 
of critical orbits can be estimated by the intersection number of special Schubert classes. 
Relations to the sh representation theory (_F ) imply that this number is the dimension 
of the space of singular vectors. 

We prove also that the spectrum of the Gaudin hamiltonians is simple for generic z. 



The Gaudin model of statistical mechanics is a completely integrable quantum spin 
chain associated to the Lie algebra s/2 = s/2(C), [G]- Denote Lx. the irreducible s/2- 
module with highest weight Xj G C. The space of states of the model is the tensor 
product L = ® ■ ■ ■ ®L\^. Associate with any Lx^ a complex number Zj^ and assume 
^1, . . . , 2;„ to be pairwise distinct numbers, z = {zi, . . . , Zn)- The Gaudin hamiltonians 
T-Ci{z), . . . ,T-Cn{z) are mutually commuting linear operators on L which are defined as 
follows, 
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where fijj is the operator which acts as the Casimir element on i-th and j-th factors of 
L and as the identity on all others. 

One of the main problems in the Gaudin model is to find common eigenvectors and 
eigenvalues of the Gaudin hamiltonians, and the algebraic Bethe Ansatz is one of the 
most effective methods for solving this problem. The idea of this method is to find some 
function with values in the space of states, and to determine a certain special value of its 
argument in such a way that the value of this function is an eigenvector. The equations 
which determine these special values of the argument are called the Bethe equations, and 
the common eigenvector corresponding to a solution of the Bethe equations, is called 
the Bethe vector, |FaTj . 



It is enough to diagonalize the Gaudin hamiltonians in the subspace of singular vectors 
of a given weight. The Bethe equations associated to the space Sing^(L) of singular 
vectors of the weight Ai + ■ ■ ■ + A„ — 2k, where A; is a positive integer, have the form 

— — — = ^ , i = l,...,k. 

1=1 ~ ^' j^i ~ 

This system is symmetric with respect to permutations of unknowns ti, . . . ,tk, and any 
(orbit of) solution t° = {t\, . . . , to this system defines an eigenvector v{tP) G Sing^(L) 
of the Gaudin hamiltonians {Hji^z)}, [H] . 

A conjecture related to Bethe Ansatz says that the Bethe vectors give a basis of the 
space of states. For the SI2 Gaudin model, the conjecture was proved in |ReVj in the 
case of generic z and generic (non-resonant) weights Ai, . . . , A„, In the case of integral 
dominant weights, which is a non-generic one, results of Sec. 9 of |ReVj imply only 
that for generic z the number of Bethe vectors is at least the dimension of the space of 
singular vectors, see also Theorem 8 in |SVj . The proof of the Bethe Ansatz conjecture 
for the s/2 Gaudin model was completed in our work with A. Varchenko |SVj . where the 
bound from above for the number of the Bethe vectors by the dimension of the space of 
singular vectors was obtained. 

In jHYj, we study the function 

k n 

i=l 1=1 l<i<j<k 

If all A; are positive integers, then this function is rational, and the Bethe equations are 
exactly the equations on the critical points of this function with non-zero critical values, 

1 9$ 

The upper bound for the number of critical orbits of the function ^{t) with non-zero 
critical values by the dimension of Sing;j,(L) is one of the main results of the paper. 
The proof is difficult, see Theorems 9-11 of |SVj . As we point out in Sec. 1.4 of jSVj . 
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the critical orbits are labeled by certain two-dimensional planes in the linear space of 
complex polynomials. This observation suggests to apply the Schubert calculus to this 
problem. In the present paper we realize this approach. 

The Wronski determinant of two polynomials in one variable defines a map from the 
Grassmannian of two-dimensional planes of the linear space of complex polynomials to 
the space of monic polynomials called the Wronski map. 

The function turns out to be the generating function of the Wronski map: the 
critical orbits label planes in the preimage under this map of the polynomial 

W{x) = {x- zi)^' ... (a; - . 

In fact, this is a reformulation of a classical result going back to Heine and Stieltjes, 
Ch. 6.8 of lEij. 

To calculate the cardinality of the preimage of the Wronski map is a problem of 
enumerative algebraic geometry, and an upper bound can be easily obtained in terms 
of the intersection number of special Schubert classes. A well-known relation between 
representation theory and the Schubert calculus ([Fj) implies that the obtained upper 
bound coincides with the dimension of Sing^(L). 

In Sec. 12 we collect known facts related to the Gaudin model. In Sec. |21 we show 
that the function defining the Bethe equations is the generating function of the Wronski 
map, and in Sec. |3 we obtain an upper bound for the number of critical orbits of the 
generating function in terms of the intersection number of special Schubert classes. 

Another conjecture related to the Gaudin model says that for generic z the Gaudin 
hamiltonians have a simple spectrum. In |ReVj it was proved that Bethe vectors are 
differ by their eigenvalues for generic Ai, . . . , A„ and real z of the form zj = , where 
s » 1. If it is known that the Bethe vectors form a basis, then the simplicity of the 
spectrum follows. 

In Sec. El we deduce the simplicity of the spectrum of the Gaudin hamiltonians for 
generic z and integral dominant s/2 weights from the relation to Fuchsian differential 
equations described in |SVj . 

The argumets presented here work and give the similar results for the sip Gaudin 
model associated with the tensor product of symmetric powers of the standard slp- 
representation, see Sec. 5 of |Slj . The link between the Bethe vectors and the Wronski 
maps seems to be useful for stufy rational maps; see jH2] where the case of rational 
functions is treated. 

The author is thankful to A. Eremenko, E. Frenkel, F. Sottile, A. Varchenko, V. Za- 
kalyukin for useful discussions, to G. Fainshtein and D. Karzovnik for stimulating com- 
ments. 
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2. Gaudin's model 

2.1. The Gaudin hamiltonians. Let mi,...,m„ be nonnegative integers, M = 
{mi, . . . ,mn). Denote L^^ the irreducible s/2-niodule with highest weight rrij. The 
space of states of the Gaudin model is the tensor product 

(1) L®^ = L^, 

Associate with any L^^ a complex number Zj, and assume z^, . . . ,z^ to be pairwise 
distinct numbers, z^ = (2;°, . . . ,z^). Let e, /, h be the standard generators of s/2, 

[e,f] = h, [h,e] = 2e, [h,f] = -2f 

and 

1^ = e (g) / + / (g) e + -/i (g) /i G s/2 ® s/2 

the Casimir element. For I < i < j < n, denote flij : L®*^ — > L*^^^ the operator 
which acts as Q on i-th and j-th factors of L®*^ and as the identity on all others. The 
hamiltonians of the Gaudin model are defined as follows, 

(2) ^^(-°) = Eio%3> ^ = 1,...,- 

2.2. Subspaces of singular vectors. Write \M\ = mi + ■ ■ • + For a nonnegative 
integer A; such that \M\ — 2k>0, define Sing^, , the subspace of singular vectors of weight 
\M\ - 2k in L®^, 

(3) Singfc = Sing;.(L^*'0 = {w e L®^^ \ ew = , hw = (|M| - 2k)w } . 
Theorem 1. ( IS VI . Theorem 5 ) We have 

k + n — 2 — rrii^ — ... — — q 
n-2 



n , 
dimSing,(L«^) = Y.{-\r Yl [ 

g=0 l<ii< --<iq<n ^ 

where we assume (^) = for a < b. 



The subspace Sing;, is an invariant subspace of the Gaudin hamiltonians for any 
< k < \M\/2, and the singular vectors generate the whole of L®^^. Therefore it is 
enough to diagonalize Gaudin hamiltonians in a given Sing;,. 

2.3. Bethe equations associated to Sing;, and z^. The Bethe equations for the 
Gaudin model associated to Sing;, and z^ have the form 

n 



1=1 ' jj^i 
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Any solution t° = (tl, • • • to this system define an eigenvector t>(t°) G Sing^ of tlie 
Gaudin liamiltonians, TCj{z^)v(t'^) = jj,jv{iP), witli eigenvalues 

k 

(5) = fi,{t ) = ^ _,o^ - 2^ Zo—fo ^ J = l,---,^- 

l^j ^^^J > i=l ^3 

2.4. The master function of the Gaudin model. Consider the following function 
in variables t = (ti, . . . ,tk) and z = {zi, . . . , Zn), 

k n 

(6) nt,^)= n (^.-^.r"^^/'nn(^^-^')"" n (^^-^^■) 



l<i<j<n i=l 1=1 ^<i<j<k 



defined on 

(7) C = C{t; z)= { t G C^ 2 G C" I 

U Zj, U ti, Zj Zg, 1 < i I < k, 1 < J ^ q < n} . 

Write 

S{t,z) = ln^(t,z) = 

k n 

^^ln(^^-2:j) - J] J]mzln(t,-2;;) + ^ 2\n(ti-tj). 

^<i<j<n 1=1 1=1 l<i<j<k 

The both functions S{t, z) and z) have clearly the same critical set in C{t] z). Denote 

Z = {z eC'lzj ^ Zg, I < 3 ^ q <n} . 

Let 71 : C —>■ Z, nit, z) = z, be the natural projection. For fixed z^ G Z, the equations 
form a critical point system of the function S{t,z^) considered as a function in t on 

n-\z'), 

dS 

(t,/)=0, z = l,...,/c, 

Oti 

and the equalities (jSj) can be re-written in the form 

yUj = -Q^{t ,z ) , J = l,...,n, 

where = is a solution to the Bethe equations. In terms of the function z), 

the Bethe equations (jH) are 

1 .|(i,.0)=0,» = l,...,i, 



■(t 

and 



1 -.0 0, 
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2.5. The Shapovalov form. Let m is a nonnegative integer. For the s/2-niodule 
with highest weight m, fix the highest weight singular vector 

Denote Bm the unique bilinear symmetric form on Lm such that 

Bra{vm, Vm) = 1, Bm{hx, y) = Bm{x, hy) , Bm{ex, y) = Bm{x, fy) 

for all x,y E Lm- The vectors fm, /fm, • • • , f^Vm are orthogonal with respect to Bm and 
form a basis of Lm- 

The bilinear symmetric form on L®^ given by 

(8) B = Bm,®---(^ Bm„ 

is called the Shapovalov form. 

Let ji, . . . , jra be integers such that < jj < mj for any 1 < i < n. Write J = 
(ji, . . .,jn) and \J\= h jn • Denote 

(9) /V/ = f ® ■ ■ ■ ® f • 

The vectors {/"^vm} are orthogonal with respect to the Shapovalov form B and provide 
a basis of the space L'®^^ . We have 

hifvM) = (|M| -2|J|)/Vf , e(/^t;M) = 0, 

i.e. the vector f'^VM is a singular vector of weight \M\ — 2\J\. The space Sing;, is 
generated by the vectors f'^VM with \J\ = k. 

2.6. Bethe vectors. For J = {ji, . . . ,jn) with integer coordinates such that < jj < 
and \J\ = k and for (t, z) in C given by ((Tj) set 



z) = . , V . Sym^ 



nn 



Jl- ■ ■ ■ 3n- \_l=l i=l ^ 

where 

SymiF(t) = ^ (t<x{i), • • • , ta{k)) 
is the sum over all permutations of ti, . . . , tfc. 

Theorem 2. ( |Re Vj ) (i) Ift^^^ is a nondegenerate critical point of the function S{t, 
then the vector 

(10) ^(t«,z°)= 5^ Aj{t^'\z')fvM 

J:\J\=k 

is an eigenvector of the operators TCi{z^), . . . , Tin^z^)- 

(a) For generic z^ , the eigenvectors v{t^'^\z^) corresponding to all critical point of the 
function S{t, generate Sing^ . < 
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The words "generic z^" mean that does not belong to a suitable proper algebraic 



The set of critical points of the function S{t, z^) is invariant with respect to the 
permutations of ti, . . . , t^, and critical points belonging to the same orbit clearly define 
the same vector. Theorem |21 gives a lower bound for the number of critical orbits of 
the function S{t, z^) by the dimension of Sing^, . In jSVj the upper bound by the same 
number was obtained. 

Theorem 3. |SVj For fixed generic , all critical points of the function S{t, z^) are 
nondegenerate, and the number of orbits of critical points is at most the dimension of 
Singfc . < 

Corollary 1. |SVj For generic z^ , the Bethe vectors of the s/2 Gaudin model form a 
basis in the space of singular vectors of a given weight. < 

The statement that all critical points of the function ^(t, z^) are nondegenerate is easy, 
see Theorem 6 in |SVj . The difficult part of |SVj was to estimate from above the number 
of critical orbits (see Theorems 9-11 in |SVj ). As we pointed out in Sec. 1.4 of |SVj . the 
orbits of critical points of the function S{t, z^) are labeled by certain two-dimensional 
planes in the linear space of complex polynomials. This observation suggests to apply 
the Schubert calculus in order to prove the following 

Proposition. The number of orbits of critical points of the function S{t, z^) is at most 
the dimension of Sing;. . 

As we will see below, within the framework of the Schubert calculus, the Proposition 
is a direct and immediate corollary of easy Theorem [7| of Sec. 14.21 



3.1. The Wronski map. Let Poly^ be the vector space of complex polynomials of 
degree at most d in one variable. Denote (j2(Poly^) the Grassmannian of two-dimensional 
planes in Poly^. The complex dimension of G2{PoIyj) is 2d — 2. 

For any V G G2(Poly^) define the degree of V as the maximal degree of its polynomials 
and the order oi V as the minimal degree of its non-zero polynomials. Let V G G2(Poly^) 
be a plane of order a and of degree b. Clearly < a < b < d. Choose in V two monic 
polynomials, F{x) and G{x), of degrees a and b respectively. They form a basis of V. 
The Wronskian of V is defined as the monic polynomial 



set of C". 



3. The generating function of the Wronski map 



Wv{x) 



F' 



{x)G{x) - F{x)G\x) 



a — b 



The following lemma is evident. 

Lemma 1. (i) The degree of Wv{x) isa + b — l<2d — 2. 

(ii) The polynomial Wv{x) does not depend on the choice of a monic basis. 

(Hi) All polynomials of degree a in V are proportional. 



< 
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Thus the mapping sending V G G'2(Poly^) to Wv{x) is a well-defined map from the 
Grassmannian G2(Poly^) to CP^'^"^. We call it the Wronski map. This is a mapping be- 
tween smooth complex algebraic varieties of the same dimension, and hence the preimage 
of any polynomial consists of a finite number of planes. On Wronski maps see 

3.2. Planes with a given Wronskian. 

Lemma 2. Any element of G2{Po\yj) with a given Wronskian is uniquely determined 
by any of its polynomial. 

Proof: Let W{x) be the Wronskian of a plane V G G'2(Poly^) and f{x) G V. Take 
any polynomial g{x) G V linearly independent with f{x). The plane V is the solution 
space of the following second order linear differential equation with respect to unknown 
function u(x), 

u{x) f{x) g{x) 
u'{x) f'{x) g'{x) =0. 
u"{x) fix) g"{x) 

The Wronskian of the polynomials f{x) and g{x) is proportional to W{x), therefore this 
equation can be re-written in the form 

W{x)u"{x) - W'{x)u'{x) + h{x)u{x) = , 

where 

-W{x)r{x) + W'{x)f'{x) 
h{x) = 77-^ ' 

as f{x) is clearly a solution to this equation. < 

We call a plane V G G2{Polj^) generic if for any xq E C there is a polynomial 
P{x) G V such that P(xo) 7^ 0. In a generic plane, the polynomials of any basis do 
not have common roots, and almost all polynomials of the bigger degree do not have 
multiple roots. 

The following lemma is evident. 

Lemma 3. Let V G G2{Polj^) be a generic plane. 

(i) If P{x) = (x — Xi) . . . {x — xi) E V is a polynomial without multiple roots, then 
Wv{xi) j^O for alll<i< I. 

(a) If Xq is a root of multiplicity /i > 1 of a polynomial Q{x) G V, then Xq is a root of 
Wv{x) of multiplicity fi — 1. < 

A generic plane V is nondegenerate if the polynomials of the smaller degree in V do 
not have multiple roots. 

Lemma 4. Let V G G2{Po\j^) be a nondegenerate plane. If the Wronskian Wv{x) has 
the form 

Wv{x) = x^'Wix) , W{0)j^O, 
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then there exists a polynomial Fq{x) & V of the form 

Fo{x) =x"'+^F{x), F(0)^0. 

Proof: Let G{x) e F be a polynomial of the smaller degree. We have G(0) 7^ 0, due 
to Lemma El Let F{x) G be a polynomial of the bigger degree. The polynomial 

F,{x) = F{x)-^^G{x) eV 

satisfies -Fo(O) = and therefore has the form 

Fo(x) = x'F(x) , F(0)^0, 

for some integer / > 1. The polynomials G and Fo form a basis of V, therefore the 
polynomial F^{x)G{x) — Fo{x)G'{x) is proportional to Wv{x). The smallest degree 
term in this polynomial is I ■ ai ■ G{0)x'-~^ , where ai is the coefficient of in Fq{x). 
Therefore / = m + 1. < 

3.3. Nondegenerate planes in Poly^ with a given Wronskian. Recall that the 
resultant Res(P, Q) of polynomials P{x) and Q{x) is an irreducible integral polynomial 
in the coefficients of P{x) and Q{x) which vanishes whenever P{x) and Q{x) have a 
common root, and the discriminant A{P) of P{x) is an irreducible integer polynomial 
in the coefficients of P{x) which vanishes whenever P{x) has a multiple root. 

Our aim is to describe the nondegenerate planes of a given order in the preimage of 
the polynomial 

(11) w{x) = {x- zir\..ix- z^r- . 

under the Wronski map. 

If plane V of order k and of degree d has the Wronskian W{x), then d = |M| + 1— A; > k. 
Let F{x) be an unknown polynomial of degree 1 < A; < \M\/2. Consider the following 
function, which appeared in unpublished notes of V. Zakalyukin related to the Wronski 
map. 

If the polynomial F{x) belongs to a plane in Poly|jvf|+i_fc with the Wronskian W{x), 
then W{x) = F'{x)Q{x) — F{x)Q'{x) for some polynomial Q{x). Differentiating gives 
W'{x) = F"{x)Q{x) - F{x)Q"{x) , and we have 

ir _ F"Q - FQ" 
W ~ F'Q - FQ' ' 

If the plane spanned by F{x) and G{x) is non-degenerate, then polynomials F and Q do 
not have common roots, polynomial F does not have multiple roots, and polynomials 



10 



I. Scherbak 



W and F do not have common roots, by Lemmas El and 01 Therefore at each root ti of 
F we get 

3.4. Theorem of Heine— Stieltjes. In order to re-write the function $ as a function 
in unknown roots of F{x), recall that if F{x) = {x — ti) . . . {x — tjS) and if W{x) is as in 
(HD), then 

k n 

A(F)= n i^^-hf^ Kes{W,F) = \['\[{U-z,r . 

We have 

k n 

(13) ^=^t-z,M)=\['\[{u-zi)-^^ n ^^^-h?- 

i=l 1=1 l<i<j<k 

The critical points of this function were studied by Heine and Stieltjes in connec- 
tion with second order linear differential equations having polynomial coefficients and 
a polynomial solution of a prescribed degree. The result of Heine and Stieltjes can be 
formulated as follows. 

Theorem 4. (Heine-Stieltjes, cf. jSzj, Ch. 6.8) Lett^ he a critical point with non-zero 
critical value of the function $(i:) given by (fT^ . Then F{x) = (a; — t?) ... (x — t°) is 
a polynomial of the smaller degree in a nondegenerate plane with the Wronskian W{x) 
given by (fTTj). 

Conversely, if F{x) = {x — t^) . . . {x — ti) is a polynomial of the smaller degree in a 
non-degenerate plane V with the Wronskian Wv{x) = W{x), then t'^ = {t^, . . . ,t1) is a 
critical point with non-zero critical value of the function $(t). < 

The function $(t) is symmetric with respect to permutations of ti,...,tk, critical 
points belonging to one orbit define the same polynomial F{x) and hence, according to 
Lemma 121 the same plane. 

3.5. Bethe vectors and nondegenerate planes. The function $(t) given by (jl3j) 
and the function \E'()f:, z) given by (jHlare differ in a factor depending only on z and M, 

^(t, ^) = n - ^,)"^»™^/2 . z, M) . 

l<i<j<n 

An easy calculation shows that ()12|) is exactly the system defining the critical points with 
non-zero critical values of the function $, considered as a function in t. The theorem of 
Heine-Stieltjes implies the following result. 

Corollary 2. There is a one-to-one correspondence between the critical orbits with non- 
zero critical value of the function (fT^ and the nondegenerate planes of order k and of 
degree \M\ -\- 1 — k having Wronskian (fTT|). < 
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Corollary 3. For fixed z & Z , M = (mi, . . . ,m„) and k, there is a one-to-one corre- 
spondence between the Bethe vectors of the s/2 Gaudin model associated to z and Sing^ 
given by (0j and nondegenerate planes of order k with the Wronskian Ull]) . < 

Corollary 4. The Bethe equations associated to z and coincide with the critical 
point system of the generating function associated to order k and the Wronskian (11111 . < 

The function $ is called the generating function of the Wronski map: for any given 
monic polynomial W{x) and any given order k < degW/2, the critical orbits of the 
function $ determine the non-degenerate planes of order k in the preimage of W{x). We 
will show below that for generic z the critical orbits determine all planes of order k in 
the preimage, see Corollary El in Sec. 14.21 

4. The preimage of a given Wronskian 

The number of nondegenerate planes with a given Wronskian can be estimated from 
above by the intersection number of Schubert classes. 

4.1. Schubert calculus (IHrHI, [E])- Let G2{d+1) = G2(C'^+^) be the Grassmannian 
variety of two-dimensional subspaces V C C^"^^. A chosen basis ei, . . . ,ed+i of C'^+^ 
defines the flag of linear subspaces 

E. : E, C E2 C ... C E^C = , 

where Ei = Spanjei, . . . , Cj} , dim Ei = i. For any integers ai and 02 such that < 02 < 
ai < d — 1 , the Schubert variety Qai,a2{E,) C G2{d -\- 1) is defined as follows, 

fiai,a2 = ^aua^iE.) = {V & G2{d + 1) | dim {V H E^-a,) > 1 , dim (\/ H Ed+l-a^) > 2 } . 

The variety ^ai,a2 = ^ai,a2{E,) is an irreducible closed subvariety of G2{d -\- 1) of the 
complex codimension ai -|- 02- 

The homology classes [^^01,02] of Schubert varieties flai,a2 are independent of the choice 
of flag, and form a basis for the integral homology of G2{d -\- 1). Denote 0"ai,a2 the 
cohomology class in H'^^°-^~^"'^\G2{d -\- 1)) whose cap product with the fundamental class 
of G2{d-\- 1) is the homology class [fia^_a2]- The classes <Jai,a2 a-^e called Schubert classes. 
They give a basis over Z for the cohomology ring of the Grassmannian. The product or 
intersection of any two Schubert classes <Jai,a2 and <Jbi,b2 has the form 

0"ai,a2 ■ <^bi,b2 = C(fll , , Cl , 02)0"^ ,C2 ; 

ci+C2=ai +02+61+62 

where C(ai, 02; &i, &2; Ci, C2) are nonnegative integers called the Littlewood-Richardson 
coefficients. 

If the sum of the codimensions of classes equals dimG2('^+l) = 2d — 2, then 
their intersection is an integer (identifying the generator of the top cohomology group 
Cd-i^-i G H'^'^^^{G2{d + 1)) with 1 G Z) called the intersection number. 



12 I. Scherbak 

When (01,02) = {q,0), < q < d — 1, the Schubert varieties flq^ are called special 
and the corresponding cohomology classes ag = cTg are called special Schubert classes. 

The Littlewood-Richardson coefficients appear as well in the decomposition of the 
tensor product of two irreducible finite-dimensional s/2-niodules into the direct sum of 
irreducible s/2-niodules. This leads to the following claim connecting Schubert calculus 
and representation theory. 

Denote Lg the irreducible sZ2-niodule with highest weight q. 

Theorem 5. jFj Let qi, . . . , qn+i be integers such that < qi < d—1 for alll < i < n+1 
and qi + ■ ■ ■ + qn+i = 2d — 2. The intersection number of the special Schubert classes, 
■ ... ■ o-q„+i , coincides with the multiplicity of the trivial sl2-module Lq in the tensor 
product Lq^^ ■ ■ ■ ® Lq^^^ . < 

Theorems ^ and El imply the following explicit formula. 

Theorem 6. Let qi, . . . , qn+i be integers such that < qi < d — 1 for all 1 < i < n + 1 
and qi + ■ ■ ■ + = 2d — 2. Then 

^q^, + ■■■ + qi^+l-d-l 

1=1 l<ii<---<i;<ra ^ n 2 

is the intersection number a „, ■ ... ■ (Ja . < 

yi yn+i 

We did not find this formula in the literature on the Schubert calculus. 

4.2. Planes with a given Wronskian and Schubert classes. Applying the Schubert 
calculus to the Wronski map, we arrive at the following result. 

Theorem 7. Let mi , . . . , m„ be positive integers, \M\ = mi + ■ ■ ■ + m„. For generic 
z & Z and for any integer k such that 1 < k < \M\ + l — k, the preimage of the polynomial 
(llip under the Wronski map consists of at most 

nondegenerate planes of order k and of degree < \M\ + 1 — k. 

Proof: Any plane of order k with the Wronskian of degree \M\ lies in G2(Poly^J^,J^_^_l_f^), as 
Lemma^shows. Therefore it is enough to consider the Wronski map on G2(Poly|M|+i-A:)- 
For any zj, define the flag J-'zj in Po\Y\J^,^_^_l_f. , 

^o{zj) C J^i{zj) C . . .J^\M\+i-k{zj),= Poly|M|+i-fe, 
where J-'i{zj) consists of the polynomials P{x) G Poly|jv/|_,_i_^ of the form 

P{x) = ai{x - Zjf'\+^-''-' + ■■■ + ao{x - Zj)^'''^^^^'' . 

We have dim J^i{zj) =2 + 1. Lemma ^ implies that the nondegenerate planes with a 
Wronskian having at Zj a root of multiplicity mj lie in the special Schubert variety 

^mj,o{J^z,) C G'2(P0ly|Ml+l-fc) • 
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The maximal possible degree of the Wronskian Wv{x) for V G Poly[^^|_,_]^_^ is clearly 
2\M\ - 2k. Denote = \M\ - 2k, the difference between 2\M\ - 2k and \M\. If moo 
is positive, it is the multiplicity ofW{x) at infinity. 

The nondegenerate planes with a Wronskian having given multiplicity moo at infinity 
lie in the special Schubert variety fimoo ,0(^00) , where J^^o is the flag 

Polyo C Polyi C ■ ■ ■ C Polj^^j^^-^_,^ . 

We conclude that the nondegenerate planes of order k which have the Wronskian (fTT|) 
lie in the intersection of special Schubert varieties 

^mi ,0 do) • 

The dimension of G'2(Poly|j\,/|+i_^,) is exactly nii + ■ ■ ■ + m„ + moo , therefore this inter- 
section consists of a finite number of planes, and the intersection number of the special 
Schubert classes 

C"mi • ••• • CTrrin " '^\M\-2k 

provides an upper bound. < 
Theorems El and together with Corollary |2l imply the Proposition of Sec. 12.61 



Corollary 5. For generic z E Z , all planes in the preimage under the Wronski map of 
the polynomial W{x) given by (|TT|l are nondegenerate. < 

5. The simplicity of the spectrum of the Gaudin hamiltonians 

5.1. Fuchsian differential equations with only polynomial solutions and the 
Wronski map. Consider a second order Fuchsian differential equation with regular 
singular points at zi,...,Zn, n > 2, and at infinity. If the exponents at zj are and 
rrij + 1, 1 < J < n, then this equation has the form 

(14) F{x)u"{x) + G{x)u'{x) + H{x)u{x) = , 

, , G(x) — m,- 
F{x)=l[{x-z,), ttH = E 

i=i j=i 

where H{x) is a polynomial of degree not greater than n — 2. On Fuchsian equations 
see Ch. 6 of [H. 

Theorem 8. Any nondegenerate plane V G G2{Po\y^ with the Wronskian ()11|) is the 

solution space of such an equation. 

Proof: Let f{x) and g{x) be non-zero monic polynomials in V such that deg/ < degg. 
As we have seen in course of the proof of Lemma El the plane V is the solution space of 
the equation 

W{x)u"{x) - W\x)u\x) + h{x)u{x) = , 
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I. Scherbak 



where 

-W{x)f"{x) + W'{x)f'{x) 

= m 

is a polynomial proportional to the Wronskian of f'{x) and g'{x). One can easily check 
that 

W'{x) _ rrij 
W{x) ^ 3; - ;Zj ' 

Moreover, if Zj is a root of W{x) of multiplicity rrij > 1, then all coefficients of the 
equation have {x — Zj)"^^~^ as a common factor, and the equation can be reduced to the 
required form (fT^ . < 

In jFrj, a correspondence between eigenvalues of the Gaudin hamiltonians and the 
s/2-opers which determine a trivial monodromy representation 

7ri(pi\{^i,...,z„,oo})^PGL2 

was established. Notice that any nondegenerate plane with a given Wronskian defines 
such an oper, and any oper with trivial monodromy defines a plane in G2{Polji^) for 
some d. 

5.2. Bethe vectors are differ by their eigenvalues. Let t'^^^ and t^"^^ be two solutions 
to the Bethe equations (jH) associated to fixed z'^ E Z and Sing^. Denote v{t^^\z'^) 
the corresponding Bethe vectors given by (Uni) and /i^^) = {fii{t^'^) , . . . , fin{t^'^)) their 
eigenvalues given by ((S)), i = 1, 2. 

Theorem 9. // = /i^^), then = vit^^^z'^). 

Proof: Let fi{x) = {x — ti^) . . . {x — t^^^) be the corresponding polynomials, i = 1, 2. 
Any of these polynomials defines a differential equation of the form ()14j) . A classical fact 
of the theory of Fuchsian equations is that if an equation of the form (jl4|) with positive 
integer mi, . . . ,m„ has a polynomial solution without multiple roots, then all solutions 
to this equation are polynomials, see Ch. 6 of R or Sec. 3.1 of [SV^ . We have 

rj - + G{x)mx) ._ 

H^x)- ^^^^^ , ^-i,2. 

Polynomials Hi{x) and H2{x) have degree at most n — 2. We will show that these 
polynomials in fact coincide. Indeed, the substitution of zj into Hi gives 

H.{z,)=m,.j^ fliz,-z,). 



An easy calculation shows that 



E 
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n. 



We conclude that polynomials Hi{x) coincide at n points, therefore Hi{x) = H2{x). 
Thus polynomials f\{x) and f2{x) are monic polynomials of the minimal degree in the 
same solution plane. Hence fi{x) = f2{x), the solutions t^^^ and t'-^-* of the Bethe 



Corollary 6. For generic G Z , the SI2 Gaudin hamiltonians (j21) have a simple spec- 
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